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Abstract —  Using  a  geometric  approach  with  ellipses, 
results  warn  that  the  autocorrelation  presence  reduces 
the  T2chart  performance,  restricting  the  chart  ability  to 
signal  a  special  cause  that  is  acting  in  the  process.  A 
control  chart  is  one  of  the  main  techniques  of  statistical 
control  of  process  and  using  this  technique  means  auto 
correlation  absence  among  the  data  of  measured  quality 
characteristics.  Some  examples  are  presented  to  illustrate 
the  detrimental  effect  of  T2 chart  when  autocorrelation  is 
present  in  the  process.  This  article  discusses  effect  of 
correlation  and  autocorrelation  in  HotellingT2  chart 
when  there  are  two  quality  characteristics X  and  Y,  whose 
autocorrelation  and  correlation  structure  is  represented 
by  a  VAR(  1 )  model. 

Keyword —  Autocorrelation;  HotellingT2 chart;  VAR(l) 
model;  statistical  process  control. 

I.  INTRODUCTION 

The  independence  among  observations  overtime 
is  one  of  the  basic  assumptions  for  control  charts  use. 
Nevertheless,  in  some  cases,  quality  characteristic 
measures  of  neighboring  items,  according  to  the  time  they 
were  produced,  may  present  some  dependence  degree 
among  the  observations.  This  phenomenon  is  called 
autocorrelation.  According  to  Sullivan[l],  many  industrial 
operations  of  streaming  present  autocorrelation  and  one  of 
the  possible  reasons  is  the  gradual  wear  of  critical 
components  of  process.  Gani[2]states  that  independence 
hypothesis  among  observations  of  a  variable  can  be 
violated  by  the  high  production  rates  that  generate 
correlation  and  dependence  among  observations  of 
neighbor  products  according  to  the  manufacturing  instant. 

The  performance  of  traditional  control  charts  can 
be  affected  by  autocorrelation  [3]. Recent  studies  present 
alternatives  to  monitor  these  types  of  processes.  Some 


authors [4]  proposed  a  control  chart  based  on  the 
minimum  euclidean  distance  to  detect  deviations  on  the 
average  in  auto  correlated  processes.  Others [5]  considered 
the  economic  design  of  control  chart  ARMA,  used  in  auto 
correlated  processes.  Franco  et  al.  [6]  used  the  AR(1) 
model  to  describe  the  oscillatory  behavior  of  the  average 

and  showed  the  optimal  parameters  of  X  chart  using  the 
Duncan  model  .Ijee  et  al. [7]  also  used  the  AR(1)  model  to 
investigate  the  effect  of  oscillatory  behavior  of  the 

average  in  the  performance  of  X  chart  with  double 
sampling.  In  order  to  reduce  the  negative  effect  of 

autocorrelation  on  the  X  chart  performance,  Hu  et  al.  [8] 
presented  a  technique  of  systematic  sampling  called  s- 
skip. 

With  advanced  technology  and  high  production 
rates  modem  systems  generated  complex  processes  that 
are  multidimensional,  that  is,  many  quality  characteristics 
are  measured  and  controlled.  Pan[9]  describes  some  of 
these  processes  that  may  have  correlated  and  auto 
correlated  observations. 

2 

The  use  ofT*"  statistics  for  monitoring  the  mean 
vector  in  multivariate  processes  is  suggested  by  Hotelling 
[10].  Seven  decades  later,  Chen  [11,12]  showed  the 
adaptive  schemes  that  improve  the  performance  of  the 

T  chart.  According  to  Hwang  et  al. [13],  autocorrelation 
increases  the  false  alarms  rate,  while  correlation  decreases 
the  control  graph  power.  The  combined  activity  of 

2 

autocorrelation  and  correlation  in  the  performance  of  T~ 
chart  is  worthy  of  investigation. 

The  aim  of  this  article  is  to  graphically  evaluate 
the  effect  of  autocorrelation  in  two  characteristics  of 
measurable  quality  X  and  Y  when  there  is  correlation 
between  observations  of  X  and  Y.The  VAR(1)  model  was 
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adopted  to  represent  the  structure  of  correlation  and 
autocorrelation.lt  was  considered  in  evaluating  that  the 
shift  in  the  mean  is  the  most  important  in  the  whole 
process  and  that  the  mean  vector  and  the  covariance 
matrix  are  known  or  estimated  accurately. 

The  article  is  organized  as  following  way: 
section  2  describes  the  model  that  represents  quality 
characteristics  when  there  is  autocorrelation  in  the 
process;  section  3  presents  some  characteristics  of  the 
Hotellingr2chart;  the  effect  of  autocorrelation  in  bivariate 
processes  is  discussed  and  evaluated  in  section  4  and, 
finally,  it  is  concluded  about  the  work  in  section  5. 

II.  AUTOREGRESSION  MODEL  AND  CROSS - 
COVARIANCE  MATRIX 

The  classical  control  procedures  in  multivariate 
processes  consider  the  basic  assumption  that  the 
observations  follow  a  multivariate  normal  distribution  and 

are  independent,  with  mean  vector  (l0and  variance- 
covariance  matrix  Z . 

Xt=H„+e,  t  =  l,2,...,T 

(1) 

where:  Xt  represents  observations  by  a  vector  of  order  p 

x  1  (p  is  the  number  of  variables);  et  are  independent 

random  vectors  of  order  p  x  7  with  multivariate  normal 
distribution  whose  mean  is  zero  and  variance -covariance 

matrix  Z(< . 

In  many  manufacturing  processes,  the 
independence  assumption  is  violated.  Autoregression 
vectors  of  first  order  -  VAR(1)  are  used  to  model 
multivariate  processes  with  temporal  correlation  among 
observations  of  the  same  variable  and  correlation  among 
observations  of  different  quality  characteristics  [14-21] 
The  VAR(1)  model  is  represented  by: 

Xf-H  =  0(XM-n)  +  £t  (2) 

being  X,  ~  Np(n,  f)  the  vector  of  observations  of 
dimension  at  instant  t  (pis  the  number  of 

variables),  fA  is  the  mean  vector, £,  is  a  random  vector 
with  independent  observations  and  multivariate  normal 
distribution  with  zero  mean  and  covariance  matrix  E  and 
o  is  a  matrix  with  autocorrelation  parameters  of  order 

{Px  P)  ■ 

ForJarrett[3],  the  Xf  cross -correlation  matrix 

f 

has  the  following  property:  T  =  +Z  .  After  some 

algebraic  manipulation,  it  is  possible  to  obtain  the 
relationship: 


Vec  r  =  (/p2 -0®oj  Vec  X  (3) 

where  ®  is  the  product  operator  ofKronecker  and  Vec  is 
the  operator  that  transforms  a  matrix  into  a  vector  by 
stacking  its  columns. 

hi  order  to  study  the  effect  of  correlation  and 

autocorrelation  of  T~  chart,  it  was  considered  a  bivariate 
process  (p=2): 


(t>  = 

a  0 

;Z  = 

'1 

P 

0  b 

_P 

1 

From  (2)  and  (3),  it  follows  that: 

g\  =(l-a2)  oXY  =  p{\-ab) 

GxY=p(l-ab)1  a2=(l -b2)1  y 
(5) 


III.  T2  CONTROL  CHART  OF  HOTELLING 

One  of  the  most  known  control  schemes  to  detect 

deviations  in  the  mean  of  multivariate  processes  I  the  T~ 
control  chart  of  Hotelling [10].  When  the  mean  vector 

(J,0=(/701,//02)  and  the  matrix  X  are  known,  the  T~ 
monitoring  statistics  of  Hotelling  is  represented  by: 

r!=(x-(i0)'  £“'(x-i*o)  (6) 

With  the  process  in  control,  T~  —  %2p  ,if  occurs 
a  special  cause  that  induces  change  in  the  mean  vector  for 
11^12)’  T~  ~  X(p,x) with  non-centrality 
parameter  A2  =  5  F—  *5  ,  being 


8  (<5j  pm,S2  P\2  M02)  ’  see 


al.[22].  When  the  mean  vector  and  the  covariance  matrix 
are  unknown  and  must  be  estimated,  the  control  limits  are 
calculated  according  to  the  monitoring  phase  [23]. 

Some  authors  [24,  25,  26,]use  the  non-centrality 

parameter  ( A" )  as  a  displacement  measure  in  the  process 
mean  vector.  In  this  case,  the  chart  performance  is 
measured  by: 


ARL  =  {l 


p r(4.«)<tc]} 


(7) 

where:  ARL  is  the  average  number  of  samples  up  to  the 
signal;  CL  is  the  control  limit  of  the  T2  chart.  The  ARL 
measures  the  average  number  of  samples  until  the 

occurrence  of  a  false  alarm  if  =  0  .  When  the  process 
is  out  of  control,  the  chart  with  the  smallest  ARL  detects 
faster  a  process  change. 
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IV.  AUTOCORRELATION  EFFECT  IN 
BIVARIATE  PROCESSES 

To  be  used  when  the  independence 
hypothesis  among  observations  of  one  or  more  quality 

characteristics  is  not  violated,  the  Hotelling  T"  chart  was 
designed.  Excluding  the  effect  of  this  hypothesis  is  quite 
detrimental  to  the  control  chart  performance.  To  study  the 
autocorrelation  effect,  it  was  considered  the  distance  from 

X  vector  to  the  (.1  mean  vector  called 
Mahalanobisdistance  [27],  The  Mahalanobis  distance  is 
given  by: 

d!=(x-^,)'z-‘(x-^,) 

(8) 


Ratio  between  cross-covariance  matrix  (F)  and 
the  elements  of  matrices  O  and  X  is  calculated  using 
equation  (3).  In  presence  of  autocorrelation  and 
correlation,  the  Mahalanobis  distance  is  given  by: 

n2=(x-n„)  rJ(x-|i0) 

(9) 

Without  generality  loss,  considering  bivariate 
a  0 
0  b 
1  P 
P  1  _ 

and  the  vector  X  =  (x\  y)  the  distance  D2  is  equivalent 


case  where  d>  = 


x  = 


and 


,  when  ( Po  ]  0-’/^o2  ^ ) 


to: 


abx 2  +2  pa2b2xy  +a2x2  -  2pa2  xy  -  ab3  y1  +abx 2  +aby2  -2pb2xy  +b2y2  -x2  +2  pxy-y2J 
(ab  - 1)  '  {-a2b2p2  +  a2b2  +  a2 p2  -  2ab  +  b2 p2  -  p2  + 1 ) 


Equation  (10)  reveals  influence  of  a  ,b  e  p  in 
distance  D2. 

If  a=b  =  0,  namely,  0  =  0  (there  is  no 
autocorrelation),  the  distance  D~  is  reduces  to: 

D 2  =(x2-2/oxy+y2)/(l-/o2) 

(11) 

When  there  is  no  autocorrelation,  that  is,  data  are 
independent  D~  -  %2fra)  •  To  evaluate  the  autocorrelation 
effect,  it  was  used  in  this  article  the  bivariate  case  and  a 
=  0.01  ( Xl-_2.,a=0m>  l  where  D2  =  10.5966. 

When  there  is  no  displacement,  the  process  is  in 
statistical  control.  The  performance  of  a  control  chart  can 
be  assessed  in  terms  of  number  of  samples  that  the  chart 
uses  to  detect  a  shift  in  the  characteristic  to  be  monitored. 
It  is  expected  in  this  case  the  signal  given  by  the  chart  to 

be  a  false  alarm  The  value  D2  =  10.5966  amounts  to  a 
false  alarm,  on  average,  for  each  200  samples  evaluated 

2 

when  using  the  Hotelling  T  chart  [28,29]. 


hi  sections  4.1  and  4.2,  there  was  the  evaluation 
of  effect  autocorrelation  of  process  in  control  (8  =  0)  and 

process  out  of  control  (5^0),  respectively  .At  the 
graphical  evaluation  of  autocorrelation  effect,  it  was 
considered  that  the  displacement  istype 

8  —  ~  P\  [  —  5  <5 2  ~  P\2~  P02  )  >  the 

occurrence  of  a  special  cause  shifts  the  mean  vector 
(10  =  (//01  =  0;//02  =  0)to  a  new  level 

Mi  =  (/An  +  ’  Pm  +  ^iu  )  ■ 

4.1  Graphical  evaluation  of  autocorrelation  effect  with 
process  in  control 

In  a  process  free  of  autocorrelation,  a=b  —  0 
and  p  =  0.7,  it  follows  that 

D2  =1.9608x2 -2.745 lxy +  1.9608y2 .  The  ellipse 

representing  the  distribution  contour  line  for  D~  = 
10.5966  is  illustrated  in  Figure  1.  Figure  2  illustrates  the 
contour  line  for  a  —b  —  0.5  and  p  =  0.7. 
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Fig.l :  Ellipse:  a  =  b  =  0  e  p  =  0.7 .  Figure 2.  Ellipse:  a  =  b  =  0.5  e  p  =  0.7. 
Source:  The  authors.  Source:  The  authors. 


For  a  and  b  S  {0.0;0.1;0.2;0.3;0.4;0.5;0.6;0.7;0.8;0.9},  it  can  be  observed  in  Figure  3  a  graphic 


demonstration  in  which  the  greater  the  autocorrelation,  the  greater  the  elliptical  area,  that  is,  if  the  process  is  in  control  and 
there  is  autocorrelation  in  variables,  it  is  necessary  to  adjust  the  graphical  control  limit  otherwise,  many  false  alarms  will 
occur. 


Source:  The  authors. 

Ellipses  in  Figure  3  represent  all  equidistant 
points,  in  Mahalanobis  distance,  from  the  origin,  if  the 
data  are  normally  distributed,.  This  suggests  that  all  these 
points  are  equally  likely  to  be  governed  by  a  multivariate 

normal  distribution  centered  at  (0,0),  for  (10=  0.  In  T~ 

Hotelling  chart,  the  control  limit  (CL)  equal  to  D~  = 
10.5966,  generates,  generally,  a  false  alarm  every  200 
samples  collected  whenfl=/?  =  0.  Not  so  when 
a -b^  0.0,  that  is,  the  average  false  alarm  rate  does 
not  correspond  to  an  alarm  every  200  samples  collected, 
even  if  it  is  used  as  the  CL  value  10.5966.  This  means  in 

practice  that,  when  one  uses  the  Hotelling  1  chart,  to 
consider  the  chart  CL  with  chi-square  distribution  with  p 


degrees  of  independence  (Z(P) )  in  autocorrelation 

presence  will  provide  with  a  false  alarms  rate  different 
than  the  desired. 

4.2  Graphical  evaluation  of  autocorrelation  effect  with 
the  process  out  of  control 

hi  Figure  4,  it  is  illustrated  a  zero-autocorrelation 
process  with  a—b  —  0  and  p  =  0.7.  The  dashed  ellipse 
with  center  (0,0)  is  a  process  in  control  and  its  equation  is 

D2  =  1.9608a-2 -2.7451xy  +  1.9608/  =10.5966 

.  The  other  ellipses  represent  the  occurrence  of  a  special 
cause  that  moves  the  mean  vector  towards  a  new  baseline: 
a)  Displacement  1  ->  =  (//01  +l;//02  +l)  ;  and: 

D 2  =  1, 960 8x2  -  2.745 lxy  - 1 . 1765 jc  + 1 .960 8 y 
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Displacement  2  ->  ^  =  (jum  +  2;  //0,  +  2)  ;  uma  causa  especial  que  desloca  o  vetor  de  medias  para 

urn  novo  patamar:  The  other  ellipses  represent  the 
occurrence  of  a  special  cause  that  moves  the  mean  vector 

Z)2  =1.9608x2  -2.745 lxy-2.3529x  +  1.9608tH2War3s^^Jefdi-7059  =  10.5966 

a)  Displacement  1  Jttj  =  ( jUQl  + 1;  jUQ1  + 1)  ;  and: 

D-  =  x2  - 1 Axy - 0. 1 8x +y2- 0. 1 8y  +  0.054 

and: 


b)  Displacement  3  ->  =(jU0l  +3;  jU{p  +3) 


D 2  =  1 .960 8x2  -  2.745  Ley  - 3.529  4x  + 1 .960 8 y 

In  Figure  5,  it  is  displayd  a  process  with 
autocorrelation  with  a  —b  =  0.7  and  p  =  0.7.  The 
dashed  ellipse  with  center  in  (0,0)  is  a  process  control  and 
its  equation  is:  D2  =  x2  —  1.4xy  +  y2  =  10.06  .  The 
value  of  10.06  was  used  in  order  to  make  a  fair 
comparison  that,  in  autocorrelation  presence,  maintains 
the  average  false  alarm  rate  equal  to  an  alarm  every  200 
samples.As  demais  elipses  representam  a  ocorrencia  de 


bJ^'^pta&mePf5^  =4^=^/01  +  2; p02  +  2)  ; 

and: 

D2  =  x2  -1.4xy-0.36x  +  y2  -0.36y +  0.216 

c)  Displacement  3  ->H1=(//0]  +  3,/i02  +3)  ; 

and: 

D2  =  x2  - 1  Axy  -  0.54x  +  y2  -  0.54y  +  0.486 


Fig. 4:  Ellipses:  a  =  b  =  0  and  p  =  0.7.  Figure5-  Ellipses:  a  =  b  =  0.7  and  p  =0.7. 
Source:  The  authors.  Source:  The  authors. 


One  can  note,  in  Figure  4,  that  in  processes 
without  autocorrelation  the  displacement  in  mean  vector 
caused  by  a  special  question  is  represented  by  ellipses 
which  depart  from  the  center  (0,0),  indicating  that  the 

T“chart,  in  this  case,  presents  superior  performance  in 
relation  to  the  case  where  autocorrelation  is  present.  In 
Figure  5,  the  ellipses  have  a  higher  resistance  to  remain 
close  to  the  center  (0,0)  when  occur  displacements  that 

disturb  the  mean  vector,  meaning  that  the  7’"  chat 
performance  is  lower  when  the  autocorrelation  is  present. 

4.3  Example 

It  was  simulated  two  types  of  bivariate  processes 
and  T"  chart  applied  to  control  the  variables  of  these 
processes.  In  the  first  case:  |4,()  =  (0;0)  ;  a=b  =  0.0;  =  0.7 


and  CL=  1 0.5966  (Af?Lo=200).  In  the  second  case: 

JJ.0  =  (0; 0)  ;  a=b=0.1;  =0.7  e  CZ=  10.06  (ARL„=2()0). 
Three  kinds  of  displacements  were  performed  at  the  mean 
vector:  8  =(1.0;  1.0),  8  =  (2.0;  2.0)  and 

8  =  (3.0; 3.0)  .Variables  observations  of  the  first  and 

second  processes  were  generated  with  the  models  of 
equation  (1)  and  equation  (2),  respectively.  Figures  6,  7 

and  8  show  the  results.  In  each  T  chart,  the  process  has 
displacements  in  the  mean  vector  from  the  sample  50.  The 
results  show  that  autocorrelation  decreases  the  chart 
power  to  detect  a  special  cause  which  operates  in  the 
mean  vector  of  the  process.  Similar  results  were 
illustrated  in  Figures  4  and  5. 
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T2  Hotelling:  a=b=0,0 


T2  Hotelling:  a=b=0,7 


Fig. 5  —HotellingT2  charts;  6  =  (e1)’  =1.0  ;  S2  =1.0). 
Source:  The  authors. 


Fig.6  -HotellingT2 charts;  5  =  (c),  =  2.0  ;  S2  =  2.0) . 
Source:  The  authors. 
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Fig.7  —Hotelling!2 charts;  8  =  (<5j  =  3.0  ; S2  =  3.0) . 


Source:  The  authors. 


V.  CONCLUSION 

It  has  been  evaluated,  in  this  article  ,the 
2 

autocorrelation  effect  in  a  T~ control  chart  since  it  is  one 
of  the  most  popular  tools  in  academia  and  industry.  The 

rji  2 

Mahalanobis  distance,  the  same  statistic  used  in  the  l 
chart,  was  used  to  represent  geometrically  a  process 
behavior  in  presence  and  absence  of  special  causes  that 
affect  the  average  value  of  the  monitored  variables. 

The  T~  chart  performance  is  affected  by  the 
autocorrelation  hypothesis  violation,  reducing  the  ability 
to  detect  deviations  in  the  mean  vector.  The  use  of 
ellipses  illustrated  how  the  data  of  a  process  behave  in 
presence  of  autocorrelation,  that  is,  the  masking 
displacement  effect  in  the  mean  vector  of  variables.  The 
displayed  examples  illustrated  reduction  that  occurs  in 

power  T  chart  to  detect  the  presence  of  a  special  cause 
that  shifts  in  the  process  mean  vector.  It  is  suggested,  in 
future  works,  presentation  of  statistics  or  techniques  that 
improve  control  chart  in  presence  of  autocorrelation. 
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